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ABSTRACT 


The  effect  of  injection  of  up  to  ten  percent  by  volume  air  in  damping 
the  amplitude  of  water  borne  pressure  pulses  was  studied.  Pressure  pulses 
were  generated  using  a  water  hammer  model.  The  water  hammer  model  con¬ 
sisted  of  a  length  of  pipe  with  water  flowing  through  it  and  a  quick  clos¬ 
ing  valve  located  downstream.  Measured  amounts  of  air  were  injected  into 
a  water  stream  and  the  valve  was  quickly  closed.  The  pres sure “time  history 
was  recorded  on  oscilloscope  photographs  at  two  locations  along  the  pipe. 

A  mathematical  model  was  developed  for  this  two  phase  model.  Small  amounts 
of  air  greatly  reduced  the  amplitude  of  the  pressure  pulses.  By  injection 
of  one  percent  by  volume  of  air  into  the  model  at  15.5  psia,  the  pressure 
amplitude  was  reduced  to  8.6  percent  of  the  theoretical  no  air  amplitude. 
The  effect  decreased  with  increased  pressure  of  the  system. 
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Introduct ion* 


In  recent  years  the  United  States  Navy  has  been  doing  a  great  deal 
of  research  and  development  in  methods  of  reducing  the  noise  level  of 
machinery  for  submarines  and  for  anti-submarine  ships.  The  fluid  systems 
of  these  ships  transmit  noises  via  pressure  pulses  in  the  systems*  con¬ 
tents,  A  similar  problem  concerning  fluid  borne  pressure  pulses  is  that 
of  vibrations  set  up  by  water  hammer  in  piping  systems.  The  objective  of 
this  study  was  to  quantitatively  determine  the  effect  of  air  in  damping 
water  borne  pressure  pulses. 
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2.  Theory. 

To  study  the  effect  of  air  in  damping  water  borne  pressure  pulses, 
it  was  necessary  to  generate  a  pressure  pulse  in  the  liquid.  For  this 
study,  the  water  hammer  model  was  used  to  generate  the  pressure  pulse. 

The  water  hammer  model  shown  in  Fig.  1  consists  essentially  of  a  head 
tank,  a  length  of  pipe,  and  a  quick  closing  valve.  The  water  is  allowed 
to  flow  through  the  pipe  with  a  velocity  head.  If  the  valve  is  quickly 
closed,  this  velocity  head  is  immediately  changed  to  pressure  head  which 
in  turn  sends  a  pressure  pulse  down  the  pipe.  This  pressure  wave  travels 
down  the  pipe  at  the  speed  of  sound  in  the  fluid  contained  in  the  system. 

Considering  the  system  shown  in  Fig.  1,  the  maximum  pressure  rise  by 
instantaneously  closing  the  valve  is  given  by  the  following  equation. 

Where:  H  =  Head  rise  above  the  Initial  static  no  flow  head,  ft.  of  water 

y 

=  Velocity  of  the  water  in  the  pipe  before  valve  closure,  ft/sec. 

g  =  Proportionality  constant  in  Newton’s  second  law  of  motion, 

°  (lbm-ft)/ (lbf-sec2). 

a  =  Velocity  of  the  pressure  pulse,  ft/sec. 

The  derivation  of  this  equation  as  derived  by  Stepanoff  / 1 3/  is  included 
in  Appendix  B.  This  equation  does  not  include  frictional  effects.  If 
the  valve  is  closed  instantaneously,  the  maximum  pressure  given  by  equa¬ 
tion  (1)  will  be  the  same  along  the  whole  length  of  the  pipe.  Obviously, 
a  valve  cannot  be  closed  instantaneously.  However,  this  same  maximum 
pressure  rise  at  a  given  location  A  along  the  pipe  will  be  attained  if 
the  valve  is  closed  in  less  time  than  it  takes  the  pressure  pulse  to  travel 
from  position  A  to  the  tank  and  to  reflect  back  to  position  A.  At  the 
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valve,  the  time  for  the  pressure  pulse  to  reflect  back  is  2h/a.  Where 
L  is  the  length  of  the  pipe  from  the  valve  to  the  tank.  At  the  center 
of  the  pipe,  the  time  for  the  pressure  pulse  to  reflect  back  is  L/a  . 

Next  to  the  tank  the  maximum  pressure  will  not  be  attained  unless  the 
valve  is  closed  instantaneously.  For  the  maximum  pressure  given  by 
equation  (1)  to  be  attained  at  a  given  point  located  at  a  distance  X 
from  the  valve,  the  closure  time  must  be  less  than  2.(L-X)/d  * 

The  speed  of  the  pressure  wave  in  a  pure  water  pipe  system  is  given 
by 

^  1 

(2) 


a  = 


Vf*  (  '/kw  4  egr) 

Where:  =  Water  density,  slugs/ft^ 

2 

K  =  Water  bulk  modulus,  Ibf/ft  . 
w 

D  =  Inside  pipe  diameter,  in. 
e  =  Pipe  wall  thickness,  in. 

E  =  Pipe  modulus  of  elasticity. 

This  equation  has  been  derived  in  Appendix  B.  This  equation  can  be  ex¬ 
tended  to  a  two  phase  air-water  system.  The  speed  of  a  pressure  pulse 
based  upon  an  isothermal  compression  of  the  air  by  the  pressure  pulse  is 
given  by 


ai  “  qj  f  +  X(  '-*)  \p  .  +e% 

q  jl  -J 


(3) 


Where:  x  =  Volume  fraction  air  at  water  inlet  conditions, 

a^.  =  Velocity  of  the  pressure  pulse  in  the  air-water  mixture 
with  isothermal  compression  of  the  air  by  the  pressure 
pulse,  ft/sec. 

a  =  Velocity  of  sound  in  air  at  the  temperature  of  the  fluid, 
a  ft/sec. 
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a  =  Velocity  of  sound  in  water  at  the  temperature  of  the  fluid, 


w 


ft/sec. 


k  =  Ratio  of  specific  heats  at  constant  pressure  to  that  at 
constant  volume  for  air, 

2 

P  =  Absolute  pressure  of  the  system,  Ibf/ft  . 

ABS 


J^a  =  Air  density,  slugs/ftP 

=  Compressibility  of  water,  ft^/lbf. 


The  speed  of  the  pressure  pulse  based  upon  an  adiabatic  compression  of  the 
air  by  the  pressure  pulse  is  given  by: 


X 

Q 


T#  +  X' 
oF 


(3a) 


Equations  (3)  and  (3a)  have  been  derived  in  Appendix  B.  It  is  noted  that 
equations  (3)  and  (3a)  reduce  to  equation  (2)  when  no  air  is  present. 
Substituting  equation  (2)  into  equation  (1)  and  multiplying  by  the  specific 
weight  of  water  (  Y^)  ,  gives  the  maximum  pressure  rise  upon  valve  closure 
with  no  air  in  the  pipe. 


HyX, 


V0 


X 


LuL. 


x/psv  (/K 


w 


(4) 


The  maximum  pressure  rise  for  an  air-water  mixture  based  upon  isothermal 
compression  of  air  is  given  by: 


P  - 


(5) 


The  term  a^  is  given  by  equation  (3).  The  derivation  of  equation  (5)  is 
in  Appendix  B.  The  same  equation  based  upon  adiabatic  compression  of  the 
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air  is  given  by: 


P  = 


Vo  Tc  a 


c 


9* 


(5a) 


The  term  a  is  given  by  equation  (3a).  The  derivation  of  equation  (5a) 

A 

is  in  Appendix  B.  The  ratio  of  equation  (5)  to  equation  (4)  gives: 


The  ratio  of  equation  (5a)  to  equation  (4)  gives: 


If  equation  (6a)  is  used,  adiabatic  conditions  are  assumed  or  if  equation 
(6)  is  used,  isothermal  conditions  are  assumed.  Normally,  sound  traveling 
in  the  atmosphere  is  adiabatic.  Karplus  / 4/  conducted  experiments  with 
air  bubbles  approximately  0.01  centimeters  in  diameter  in  an  air-water 
mixture  and  found  that  his  experimental  data  was  best  correlated  by  assum¬ 
ing  an  isothermal  compression.  The  method  used  to  inject  air  for  this 
investigation  did  not  give  the  small  bubble  distribution  Karplus  obtained. 
In  Karplus1  method,  there  was  a  high  heat  transfer  surface  per  unit  bubble 
volume.  Therefore,  the  compression  was  nearly  isothermal  in  the  water. 
With  larger  bubbles  there  is  less  surface  per  unit  bubble  volume  and  the 
situation  is  probably  best  described  by  the  adiabatic  equation  (6a). 

Equation  (6a)  provides  the  fundamental  theory  for  the  damping  of 
pressure  pulses  in  water  by  air  injection.  By  measuring  the  air  injection 
rate,  the  water  flow  rate,  and  the  pressure-time  history  at  a  given  loca¬ 
tion  in  the  system  of  Fig.  1,  the  effect  of  air  in  damping  a  pressure 
pulse  in  water  is  obtained.  In  this  report,  the  experimental  pressures 
obtained  are  compared  with  the  theoretical  relationship  of  equation  (6a). 
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3.  Equipment. 

To  study  the  effect  of  air  in  damping  water  borne  pressure  pulses, 
and  to  verify  the  applicability  of  equation  (6a)  to  the  water  hammer 
model,  the  system  shown  in  Fig.  1  was  constructed. 

The  system  consisted  of  1)  a  head  tank,  a  pipe  and  a  weigh  tank 
system  2)  an  air  injection  system  3)  a  valve  and  closure  system  and 
4)  a  pressure  recording  system.  A  general  description  for  each  part  of  the 
system  will  be  presented  here.  Exact  detailed  equipment  specifications  are 
included  in  Appendix  C. 

The  head  tank  was  a  tank  two  feet  in  diameter  with  four  constant  head 
overflow  sections.  These  overflow  sections  were  located  nominally  at  two, 
four,  six  and  eight  feet  above  the  main  tank  outlet  line.  The  main  pipe 
length  consisted  of  59  feet  6  and  7/8  inches  of  two  inch  extra  strong 
black  steel  pipe.  This  pipe  was  supported  at  ten  feet  intervals  by  a  timber 
structure.  At  the  main  pipe  discharge,  a  20  inch  by  28  inch  by  20  inch 
weigh  tank,  platform  scale  and  electric  timer  were  used  to  measure  the 
weight  flow  rate  of  the  water.  Water  for  the  system  was  supplied  from  the 
city  water  supply. 

The  air  injection  system  contained  a  Fischer  and  Porter  Company  roto- 
meter  with  a  glass  ball  float  for  air  flow  measurement.  This  rotometer 
had  a  maximum  capacity  of  0.325  cubic  feet  per  minute  at  14.7  psia  and 
70°F.  Compressed  air  from  a  nominal  100  psig  compressor  was  reduced  by  a 
regulating  valve  to  approximately  fifteen  psig,  A  mercury  manometer  was 
used  to  measure  the  air  pressure,  and  the  air  temperature  was  measured 
with  a  copper  constantan  thermocouple  junction  upstream  of  the  rotometer. 
The  air  was  transferred  from  the  regulating  valve  to  the  point  of  injection 
through  one  quarter  inch  vinyl  tubing.  The  air  was  regulated  by  a  manual 
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control  valve  located  immediately  following  the  rotometer.  Fig.  1  shows 
a  diagram  of  the  system.  Fig.  2  shows  a  photograph  of  the  component  parts 
of  the  air  injection  system.  The  air  was  injected  into  the  pipe  near  the 
tank  through  a  system  of  three  circles  of  small  holes  drilled  equally 
spaced  around  the  circumference  of  the  pipe.  Details  are  included  in 
Appendix  C.  A  concentric  plexiglass  cylinder  with  "0"  ring  seals  surrounds 
these  holes.  The  system  was  designed  so  that  any  one,  any  two,  or  any  three 
of  the  circles  of  holes  could  be  used  to  inject  air  at  any  one  time.  Fig. 

5  shows  a  cross  sectional  view  of  the  air  inject  annular  ring.  This  method 
for  air  injection  is  very  simple  and  could  easily  be  adapted  in  practical 
applications . 

If  the  same  maximum  pressure  (neglecting  friction  losses)  was  to  be 
obtained  along  one  half  the  length  of  the  pipe,  it  was  necessary  to  close 
the  valve  in  less  than  one  half  the  pipe  period  This  was  obtained 

by  the  use  of  a  modified  two-inch  Crane  quick-closing  valve.  A  steel  cy¬ 
linder  two  inches  in  diameter  and  one  inch  high  was  screwed  onto  the  top 
of  the  stem.  The  valve  was  closed  by  dropping  a  seven  pound  weight  from 
nine  feet  onto  the  cylinder.  A  spring  loaded  catch  was  provided  to  keep 
the  valve  closed  when  the  valve  was  in  the  shut  position.  An  isolated 
stop  stand  was  provided  for  stopping  the  valve.  This  step  stand  minimized 
the  effect  of  the  valve  stopping  impact  from  being  transmitted  to  the  pipe 
and  to  the  valve  body.  A  quick  release  mechanism  was  provided  to  release 
the  weight  which  was  guided  down  through  a  perforated  aluminum  tube,  one 
and  one-quarter  inches  in  diameter.  Fig.  3  shows  a  photograph  of  the 
valve  and  closure  system.  In  Appendix  G,  calculations  are  provided  which 
were  used  to  estimate  the  closure  time  of  the  valve.  The  closure  time  of 
the  valve  was  measured  with  an  oscilloscope  and  a  Polaroid  camera.  The 
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time  sweep  was  triggered  with  a  photoelectric  cell.  The  closure  was 
marked  by  a  45  volt  direct  current  signal  on  the  vertical  plates  of  the 
oscilloscope.  The  distance  that  the  sweep  moved  on  the  oscilloscope 
from  trigger  to  closure,  was  photographed  with  3000  speed  Polaroid  film. 

The  closure  time  for  the  seven  pound  weight  from  nine  feet  was  0,007  +  0,003 
seconds.  This  agrees  closely  with  the  values  calculated  in  Appendix  G. 

A  Kistler  model  401  quartz  pressure  transducer  was  installed  at  three 
locations  along  the  pipe,  A  Kistler  407-5  low  noise  cable  was  used  to 
connect  the  transducer  to  the  Kistler  model  651  piezo-calibrator  amplifier, 

A  shielded  cable  was  used  to  connect  the  output  of  the  piezo-calibrator 
amplifier  to  a  Hewlett  Packard  model  130  oscilloscope.  An  oscilloscope 
Polaroid  camera  was  used  with  Polaroid  3000  speed  film  to  take  a  time 
exposure  of  the  pressure- time  history  for  each  run.  When  the  valve  was 
completely  closed,  a  45  volt  direct  current  circuit  was  completed.  This 
triggered  the  oscilloscope  time  sweep.  Fig.  4  is  a  photograph  of  the  com¬ 
ponent  parts  of  the  pressure  recording  system. 


8 


4. 


Procedure . 


The  objective  of  the  experimental  part  of  this  work  was  to  verify 
the  applicability  of  the  theoretical  model. 

The  pressure-time  history  was  recorded  at  three  locations  along  the 
pipe.  Fig.  13  shows  these  locations.  They  were  numbered  starting  from 
the  valve,  one,  two,  and  three  respectively.  The  pressure-time  history 
was  recorded  at  each  transducer  location  for  five  different  air-water 
mixtures  at  each  nominal  head  in  the  head  tank  of  two*  four,  six  and 
eight  feet.  Since  the  effect  of  small  amounts  of  air  was  of  interest, 
air  rates  were  chosen  such  that  the  volume  percent  air  in  the  system  was 
always  less  than  ten  percent. 

Because  of  the  effect  of  as  little  as  0.1  percent  by  volume  air  in 
the  water  for  the  no  air  runs,  it  was  difficult  to  get  good  results  with 
water  taken  directly  from  the  water  main.  For  these  runs,  it  was  found 
that  better  results  were  obtained  if  the  water  was  left  in  the  tank  over¬ 
night. 

The  transducer  was  calibrated  using  an  Ashcroft  dead-weight  tester. 

The  calibration  curves  are  given  in  Appendix  C,  Fig.  15  and  16. 

A  calibration  curve  for  the  rotometer  was  available  from  the  manu¬ 
facturer.  Four  points  on  this  curve  were  checked  with  a  positive  displace¬ 
ment  wet  gas  meter.  The  points  checked  closely,  and  the  curve  was  used. 

The  detailed  procedure  for  each  run  is  given  in  Appendix  D  and  a 
sample  data  sheet  is  shown  in  Appendix  E. 

Only  one  pressure  transducer  was  used.  This  transducer  was  moved 
to  the  desired  location  for  each  run. 
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5. 


Results . 


The  results  have  been  presented  in  the  form  of  tables,  photographs, 
and  curves.  In  Table  1  ’'transducer  locations"  are  numbered  as  shown  in 
Appendix  A,  Fig.  13.  The  "head”  refers  to  the  head  of  water  in  the  head 
tank.  The  values  given  are  nominal,  and  the  exact  values  of  the  heads  are 
shown  in  Appendix  A,  Fig.  13.  The  term  Y~|  is  the  volume  percent  of  air 
at  the  water  inlet  conditions  of  temperature  and  pressure.  The  term 
is  the  weight  percent  of  air  in  the  mixture.  P  is  the  maximum  pressure 
rise  above  the  static  pressure  with  no  flow  in  the  pipe.  The  ratio  P/P^, 
is  the  ratio  of  the  maximum  pressure  obtained  when  the  valve  is  closed  to 
the  theoretical  maximum  pressure  with  no  air,  equation  (4). 

Results  are  tabulated  for  transducer  locations  one  and  two.  Meaningful 
results  were  unattainable  from  location  three.  The  reasons  for  this  will  be 
discussed  later. 

Fig.  6  shows  a  series  of  photographs  for  nominal  heads  of  two  and  four 
feet  of  water  and  Fig.  7  shows  the  series  of  photographs  for  nominal  heads 
of  six  and  eight  feet  of  water.  Each  row  of  photographs  is  for  a  given 
head.  In  the  first  photograph  no  air  was  injected  into  the  water.  The 
amount  of  air  increases  moving  from  left  to  right.  The  camera  used  to  take 
the  photographs  takes  the  pictures  such  that  the  sweep  of  the  oscilloscope 
moves  from  right  to  left  rather  than  left  to  right  as  the  oscilloscope 
sweep  actually  moves.  The  initial  point  in  Fig.  6  and  7  has  been  marked  by 
lines.  The  horizontal  arrows  point  in  the  direction  that  the  sweep  moves 
along  the  horizontal  scale.  The  vertical  arrows  point  in  the  direction  of 
increasing  pressure.  In  the  photographs,  it  should  be  noted  that  the  pres¬ 
sure  and  time  scales  are  different  when  air  Is  Injected  from  those  when  no 
air  is  injected. 
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In  Fig.  8  through  11,  the  volume  percent  of  air  was  based  upon 
the  water  conditions  of  temperature  and  pressure  at  the  inlet  to  the  main 
run  of  pipe.  Fig.  8,  9,  10,  and  11  are  respectively  for  nominal  heads  of 
two,  four,  six  and  eight  feet  of  water.  Both  the  theoretical  adiabatic 
(Equation  (6a))  and  theoretical  isothermal  (Equation  (6))  curves  have  been 
drawn  in.  The  experimental  points  for  transducer  locations  one  and  two  are 
plotted  and  marked  with  symbols. 
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6. 


Discussion  of  Results, 


Transducer  location  three  could  not  be  used  to  obtain  meaningful 
data  because  the  valve  could  not  close  fast  enough.  The  valve  would 
have  had  to  close  in  less  than  0.0015  seconds  to  get  the  maximum  pressure 

a' 

without  interference  of  the  reflected  pulse  returning  from  the  tank.  Hence 
only  data  from  transducer  locations  one  and  two  are  presented.  This  sub¬ 
ject  was  previously  discussed  in  the  theory  section. 

As  the  tabulated  results  show,  two  runs  were  made  for  each  set  of  run 
conditions.  This  was  done  as  a  means  of  checking  the  accuracy  and  consist¬ 
ency  of  the  results. 

The  weight  percent  of  air  in  the  two  phase  water-air  mixture  is  includ 
ed  in  the  tabulated  results  to  show  the  small  weight  fraction  of  air  pre¬ 
sent  in  the  air-water  mixtures. 

It  was  impossible  to  obtain  meaningful  results  at  transducer  location 
one  with  no  air  in  the  water.  The  high  frequency  pressure  pulses  set  up 
by  the  valve  stem  and  gate  vibrating  made  it  impossible  to  distinguish 
between  the  pulse  created  by  stopping  the  water  and  that  created  by  the 
vibrating  valve  parts.  With  a  small  amount  of  air,  these  high  frequency 
pulses  were  damped  out  and  meaningful  results  were  obtained.  The  effect 
of  air  in  damping  the  high  frequency  noise  is  seen  in  Fig.  6  and  7. 

The  pressure-t ime  photographs  of  Fig. 6  and  7  show  how  the  pressure¬ 
time  history  at  transducer  location  two  appeared  on  the  oscilloscope.  The 
photographs  show  that  in  the  region  observed;  1)  the  greater  the  percent 
of  air  present  in  the  mixture,  the  smaller  the  pressure  amplitude  and  the 
slower  the  pressure  pulse  moves  down  the  pipe.  The  slowing  rate  of  the 
pressure  pulse  is  shown  by  the  increasing  time  duration  of  the  pulse  as 
the  air  percent  increases.  2)  Small  amounts  of  air  have  a  large  effect  in 
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reducing  the  amplitude  of  the  pressure  pulse  and  in  decreasing  the  pressure 
pulse  speed. 

The  photographs  taken  at  transducer  location  one  show  the  same  effects 
as  the  photographs  taken  at  transducer  location  two.  As  a  result,  the 
photographs  taken  at  transducer  location  one  are  not  included  in  this  re¬ 
port.  Fig.  6  and  7  show  the  photographs  taken  at  transducer  location  two. 
The  results  taken  from  the  photographs  are  tabulated  in  Table  1. 

In  Fig.  8  through  11,  the  theoretical  curves  include  no  frictional 
losses.  Because  there  are  no  frictional  losses  included  in  the  theoretical 
curves,  it  would  be  expected  that  the  experimental  results  should  fall 
slightly  below  the  theoretical  curve.  Thi*  is  because  the  frictional 
losses  for  a  two  phase  mixture  are  greater  f\an  those  for  a  single  phase. 
Martinelli  experimentally  showed  this.  / 1 0/  As  discussed  in  the  theory 

section,  the  compression  of  the  air  by  the  pressure  pulse  was  probably  best 
described  as  adiabatic  because  the  method  of  injection  allowed  for  large 
bubbles  which  gave  a  small  heat  transfer  area  per  unit  volume.  The  pres¬ 
sure  pulse  travels  so  fast  that  the  air  is  compressed  for  a  very  short 
time,  consequently  this  permits  very  little  time  for  heat  transfer. 

For  the  runs  when  no  air  was  injected  into  the  water,  values  of  81.7, 
71.7,  74.7  and  79.9  percent  of  the  theoretical  maximum  pressure  were  obtain¬ 
ed  at  transducer  location  two  for  nominal  heads  of  two,  four,  six,  and 
eight  feet  respectively.  Entrained  air  in  the  water  of  about  0.02  percent 
by  volume  would  cause  this  effect.  The  entrained  air  was  under  pressure 
of  about  100  psia  in  the  United  States  Naval  Postgraduate  School  system. 

When  the  water  pressure  was  reduced,  the  air  expanded  and  the  volume  per¬ 
cent  increased.  A  glass  of  water  will  have  air  bubbles  on  the  glass  sur¬ 
face  if  allowed  to  sit  overnight. 
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By  allowing  the  water  to  stand  about  fifteen  hours  in  the  tank,  the 
maximum  pressure  obtained  when  no  air  was  injected  increased  from  approxi¬ 
mately  50  percent  to  about  80  percent  of  the  theoretical  maximum.  This 
indicates  that  the  runs  in  which  no  air  was  injected  had  enough  entrained 
air  in  the  water  that  the  exact  air  content  was  uncertain.  A  factor  of 
0.02  percent  by  volume  makes  a  large  difference  in  the  maximum  pressure 
obtained  in  the  region  approaching  zero  percent  of  air.  This  same  per¬ 
cent  of  uncertainty  makes  very  little  difference  in  the  maximum  pressure 
in  the  region  greater  than  about  0.5  percent  of  air.  Therefore,  this 
entrained  air  did  not  effect  the  results  appreciably  for  values  greater 
than  0.5  percent  of  air. 

The  curves  of  Fig.  8  through  12  were  plotted  on  a  scale  in  which  the 

low  ranges  of  pressure  and  volume  percent  were  expanded.  This  was  done  to 

better  show  the  region  of  greatest  change.  A  log  scale  was  used  in  which 

for  the  ordinate  In  (  P /PNr  +  0.1)  was  plotted  and  for  the  abscissa 

In  (n  +  1  .0)  was  plotted.  This  was  done  so  that  it  would  be  possible  to 

have  zeros  for  P/P^,  anc*  ^  cn  the  scales.  The  results  indicate  that  the 

points  fall  close  to  the  theoretical  curves  in  the  regions  greater  than  two 

percent  of  air.  In  the  region  of  less  than  two  percent  of  air,  the  results 

obtained  are  more  erratic.  The  uncertainty  in  the  region  of  0.5  percent 

by  volume  air  was  plus  or  minus  about  ten  percent  in  the  volume  percent 

term,  and  about  plus  or  minus  11.0  percent  in  the  P/Pv_  term.  Even  consider 

NT 

ing  these  uncertainties,  several  of  the  points  in  the  region  less  than  two 
percent  by  volume  air  do  not  fall  on  the  theoretical  curve.  In  the  region 
of  greater  than  two  percent  by  volume  air  the  consistency  of  the  data  vas 
much  better.  At  four  percent  by  volume  air,  the  uncertainty  of  the  volume 
percent  term  (Y[  )  was  about  plus  or  minus  0.8  percent  and  that  of  the  P/P^ 
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term  was  about  plus  or  minus  17.0  percent.  These  uncertainties  are 
based  upon  the  uncertainty  analysis  of  runs  3-9-2  and  3-9-10  given  in 
Appendix  H. 

Fig.  12  shows  the  effect  of  air  in  damping  pressure  pulses  in  systems 
at  pressure  up  to  500  psia.  This  figure  shows  curves  of 
for  pressures  of  500,  400,  300,  200,  100,  and  15.5  psia.  The  volume  per¬ 
cent  of  air  at  inlet  conditions  varies  from  zero  to  ten  percent.  These 
curves  show  that  the  lower  the  pressure  in  which  the  pressure  pulse  is 
being  transmitted,  the  greater  is  the  effect  of  air  in  reducing  the  pressure 
amplitude.  As  an  example,  at  one  percent  by  volume  of  air  at  15.5  psia, 
the  pressure  would  be  reduced  to  8.6  percent  of  the  original  pressure. 

At  100,  200,  300,  400,  and  500  psi  the  pressure  amplitude  reduction  would 
be  to  20.8,  28.7,  34.4,  38.9,  and  42.5  percent  respectively  of  the  original 
pressure  amplitude. 

As  an  example  of  the  usefulness  of  air  in  reducing  pressure  pulses 
and  also  to  get  an  idea  of  the  amount  of  air  that  would  be  required,  con¬ 
sider  the  following  example. 

Example  1.  Pure  water  is  flowing  through  a  two  inch  extra  strong 
pipe  at  the  rate  of  500  gpm.  The  system  has  an  operating  pressure 
of  100  psia  and  an  operating  temperature  cf  70°F.  Fluid  borne  noises 
are  being  transmitted  in  this  system.  It  is  desired  to  reduce  the 
amplitude  of  these  pressure  pulses  to  20  percent  of  the  original 
amplitude . 

Fig.  12  shows  that  for  the  pressure  amplitude  to  be  20  percent 
of  the  original  amplitude,  a  1.05  percent  by  volume  mixture  must  be 
made  by  the  injection  of  air. 
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.500,9  prfl- 

■7.48  gal/ffi 


—  46-9  4  ft  /min  water 


Let  x  equal  the  air  required  in  ft  /min  at  100  psia  and  70°F. 
Then 

0.0105  = 


100X 


X  +  66.84 

x  =  0.7093  ft3/min  at  100  psia  and  70°F 

3 

This  equals  4.825  ft  /min  at  14.7  psia  and  70°F. 

The  weight  of  air  required  would  be  0.3628  lbm/min. 
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7 . 


Conclusions , 


Based  upon  the  experimental  work  and  development  of  equations  (6) 
and  (6a)  the  following  conclusions  have  been  made. 

1.  Small  amounts  of  air  in  a  water  system  greatly  reduce  the 
amplitude  of  a  pressure  pulse  transmitted  in  the  water.  In  a  15.5  psia 
water  system  with  one  percent  by  volume  air,  the  pressure  is  reduced  to 
8.6  percent  of  the  theoretical  no-air  amplitude. 

2.  The  effect  of  air  in  damping  pressure  pulses  decreases  as 
the  pressure  of  the  system  increases  as  demonstrated  by  Fig.  12. 

3.  There  was  a  reasonable  correlation  of  experimental  results 
with  theory  as  shown  in  Fig.  8  through  11. 
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Fig.  1  System  Diagram 


KEY  TO  SYSTEM  DIAGRAM 


1.  Water  storage  tank 

2.  Two  inch  water  fill  line 

3.  Storage  tank  overflow  piping  system 

4.  Overflow  drain  tank 

5.  Air  injection  system 

6.  Two  inch  pipe  coupling 

7.  Fischer  and  Porter  rotometer 

8.  Potentiometer  reference  ice  bath 

9.  Air  manifold 

10.  Inlet  air  mercury  manometer 

11.  Thermocouple  potentiometer 

12.  Kistler  Model  401  piezoelectric  pressure  transducer 

13.  Kistler  #470  shielded  cable 

14.  Kistler  Model  #651  piezo-calibrator 

15.  Shielded  calibrator-oscilloscope  cable 

16.  Hewlett  Packard  Model  130  oscilloscope 

17.  Polaroid  oscilloscope  camera 

18.  Falling  weight  guide  tube 

19.  Seven  pound  weight 

20.  Modified  Crane  two  inch  quick  closing  valve 

21.  Quick  closing  valve  stop 

22.  Water  discharge  mercury  thermometer 

23.  Dayton  scales 

24.  Weigh  tank 

25.  Air  inlet  themocouple 


19 


26.  Water  inlet  Weston  thermometer 

27.  Standard  ASA  two  inch  extra  strong  steel  pipe 

28.  One  quarter  inch  vinyl  tubing 

29.  Two  inch  gate  valves 

30.  Air  regulating  valve 

31.  Quick  closing  air  stop  valve 

32.  Air  flow  regulating  valve 

33.  Water  inlet  regulating  valve 

34.  Weigh  tank  two  inch  discharge  gate  valves 

35.  Oscilloscope  sweep  trigger 
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FIG.  2.  Photograph  of  Air  Injection  System 


t?io.  Fhotogranh  of  Valve  and  Closure  Svstem 


KEY  TO  EQUIPMENT  PHOTOGRAPHS 


Item 


1. 


2. 


3. 


A. 


5. 


6. 


7. 


8. 


9. 


10. 


11. 

12. 


13. 


14. 


15. 


Description 
Thermocouple  potentiometer 
Fisher  and  Porter  rotometer 
Potentiometer  reference  ice  bath 
Inlet  air  mercury  manometer 
Air  injection  system 

Hewlett  Packard  model  130  oscilloscope 
Kistler  model  #6 5 1  piezo-calibrator 
Kistler  #470  shielded  cable 

Kistler  model  401  piezo-electric  pressure  transducer 

Modified  Crane  quick  closing  valve 

Falling  weight  guide  tube  with  weight  in  it 

Weigh  tank 

Dayton  scales 

Weigh  tank  discharge  valve 
Electric  timer 
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Fig.  5  Cross  Sectional  View  of  Air  Injection  Annular  Ring 
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FIG.  6  PRESS'*!^  TIME  'HOTOGRAFHS  PAPT  1 
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Table  1.  Tabulated  Results 


Run  No. 

Trans¬ 

ducer 

Location 

Head 

ft. 

Yl 

7. 

Y 

'  WT 

7oxl0 

P 

psi 

P/P_ 

NT 

3-21-32 

1 

2 

1.09 

1.42 

10.6 

0.0498 

3-21-33 

1 

2 

1.09 

1.42 

10.6 

0.0498 

3-21-34 

1 

2 

3.05 

4.08 

10.6 

0.0498 

3-21-35 

1 

2 

3.10 

4.14 

11.2 

0.0526 

3-21-36 

1 

2 

5.14 

7.02 

8.8 

0.0413 

3-21-37 

1 

2 

5.17 

7.06 

8.2 

0.0385 

3-21-38 

1 

2 

7.19 

10.05 

7.3 

0.0343 

3-21-39 

1 

2 

7.19 

10.04 

7.3 

0.0343 

3-21-40 

1 

2 

9.23 

13.16 

6.4 

0.0300 

3-21-41 

1 

2 

9.28 

13.23 

6 . 4 

0.0300 

3-21-18 

1 

4 

0.664 

1.83 

22.4 

0.0704 

3-21-19 

1 

4 

0.664 

1.83 

22.4 

0.0704 

3-21-20 

1 

4 

2.02 

2.83 

14.6 

0.0462 

3-21-21 

1 

4 

2.01 

2.80 

15.8 

0.0497 

3-21-22 

1 

4 

3.37 

4.78 

11.1 

0.0349 

3-21-23 

1 

4 

3.37 

4.77 

11.9 

0.0374 

3-21-25 

1 

4 

4.77 

6.85 

10.4 

0.0327 

3-21-26 

1 

4 

6.38 

9.31 

10.4 

0.0327 

3-21-27 

1 

4 

6.38 

9.31 

10.4 

0.0327 

3-21-3 

1 

6 

0.538 

0.774 

30.6 

0.0765 

3-21-4 

1 

6 

0.536 

0.769 

28.8 

0.0720 

3-21-5 

1 

6 

1.53 

2.23 

26.4 

0.0660 

3-21-6 

1 

6 

1.53 

2.23 

24.6 

0.0615 
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Table  1.  Tabulated  Results  (continued) 


Run  No. 

: 

Trans¬ 

ducer 

Location 

Head 

ft. 

n 

% 

W  WT 

%xl0"3 

P 

P/Pxw. 

NT 

3-21-7 

1 

6 

2.60 

3.84 

18.0 

0.0450 

3-21-8 

1 

6 

2.60 

3.84 

18.0 

0.0450 

3-21-13 

1 

6 

3.68 

5.49 

15.6 

0.0390 

3-21-14 

1 

6 

3.68 

5.49 

15.6 

0.0390 

3-21-15 

1 

6 

4.76 

7.18 

12.0 

0.0300 

3-21-16 

1 

6 

4.80 

8.07 

14.4 

0.0360 

3-13-3 

1 

8 

0.475 

0.731 

35.8 

0.0773 

3-13-4 

1 

8 

0.475 

0.731 

41.7 

0.0901 

3-13-5 

1 

8 

1.28 

2.04 

33.5 

0.0724 

3-13-6 

1 

8 

1.28 

2.01 

34.6 

0.0747 

3-13-7 

1 

8 

2.14 

3.35 

25.7 

0.0555 

3-13-8 

1 

8 

2.14 

3.35 

29.3 

0.0633 

3-13-9 

1 

8 

3.01 

4.77 

19.1 

0.0413 

3-13-10 

1 

8 

3.01 

4.77 

20.9 

0.0451 

3-13-11 

1 

8 

3.91 

6.22 

17.9 

0.0387 

3-13-12 

1 

8 

3.91 

6.22 

19.1 

0.0413 

2-27-4 

2 

2 

0 

0 

174.0 

0.817 

4-1-4 

2 

2 

0.294 

0.382 

43.1 

0.202 

4-1-5 

2 

2 

0.294 

0.382 

45.4 

0.213 

2-28-4 

2 

2 

1.20 

1.65 

8.68 

0.0410 

2-28-5 

2 

2 

1.19 

1.58 

7.78 

0.0365 

2-28-13 

2 

2 

3.13 

4.25 

8.07 

0.0379 

2-28-14 

2 

2 

3.10 

4.22 

9.97 

0.0468 
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Table  1  Tabulated  Results  (continued) 


Run  No. 

Trans¬ 

ducer 

Location 

Head 

ft. 

n 

% 

r/ 

'  WT  - 

7»xlU 

P 

P/P^ 

NT 

2-28-15 

2 

2 

5.20 

7.23 

5.67 

0.0266 

2-28-16 

2 

2 

5.18 

7.22 

6.27 

0.0294 

2-28-17 

2 

2 

7.20 

10.2 

6.03 

0.0283 

2-28-18 

2 

2 

7.28 

10.1 

5.37 

0.0252 

2-28-19 

2 

2 

9.54 

13.9 

5.97 

0.0280 

3-27-3 

2 

4 

0 

0 

228 

0.717 

4-1-6 

2 

4 

0.277 

0.257 

105.3 

0.331 

4-1-7 

2 

4 

0.277 

0.257 

84.1 

0.261 

2-28-26 

2 

4 

0.769 

1.11 

37.9 

0.118 

2-28-28 

2 

4 

2.087 

3.12 

19.9 

0.0626 

2-28-29 

2 

4 

2.087 

3.12 

15.7 

0.0494 

2-28-30 

2 

4 

3.49 

5.18 

12,1 

0.0380 

2-28-31 

2 

4 

3.51 

5.20 

12.7 

0.0399 

2-28-32 

2 

4 

5.02 

7.55 

9.10 

0.0286 

2-28-33 

2 

4 

4.99 

7.48 

12.1 

0.0380 

2-28-34 

2 

4 

6.24 

9.51 

9.1 

0.0286 

2-28-25 

2 

4 

6.20 

9.67 

10.9 

0.0343 

3-27-2 

2 

6 

0 

0 

228.1 

0.717 

4-1-8 

2 

6 

0.136 

0,207 

161.8 

0,404 

3-3-4 

2 

6 

0.596 

0.877 

28.6 

0.0715 

3-3-5 

2 

6 

0.602 

0.893 

38,8 

0,0956 

3-3-6 

2 

6 

1.57 

2.33 

16.9 

0.0422 

3-3-7 

2 

6 

1.59 

2.41 

26.2 

0.0655 

3-3-8 

2 

6 

2.59 

3.90 

22.0 

0.0550 
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Table  1  Tabulated  Results  (continued) 


Run  No. 

Trans¬ 

ducer 

Location 

Head 

ft. 

n 

% 

3-3-9 

2 

6 

2. 

3-3-10 

2 

6 

3. 

3-3-11 

2 

6 

3. 

3-3-12 

2 

6 

4. 

3-3-13 

2 

6 

4. 

3-27-1 

2 

8 

0 

4-1-10 

2 

8 

0. 

3-9-2 

2 

8 

0. 

3-9-3 

2 

8 

0. 

3-9-5 

2 

8 

1. 

3-9-6 

2 

8 

2. 

3-9-7 

2 

8 

2. 

3-9-8 

2 

8 

3. 

3  9-9 

2 

8 

3. 

3-9-10 

2 

8 

3. 

3-9-11 

2 

8 

3. 

7x10 

P 

P/P*™ 

NT 

3.89 

16.0 

0.0400 

5.51 

17.8 

0.0445 

5.51 

17.8 

0.0445 

7.24 

15.4 

0.0385 

7.32 

10.6 

0.0265 

0 

370 

0.799 

0.179 

356.9 

0.771 

0.750 

51.9 

0.112 

0.747 

42.3 

0.0913 

2.12 

48.9 

0.106 

3.41 

26.7 

0.0576 

3.45 

27.9 

0.0602 

4.76 

19.5 

0.0421 

4.76 

21.3 

0.0460 

6.37 

17.1 

0.0462 

6.37 

15.5 

0.0462 

58 

61 

61 

70 

76 

117 

485 

484 

26 

13 

15 

01 

01 

89 

89 
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APPENDIX  A 
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Fig.  13  Pipe  and  Tank  Dimensions. 


APPENDIX  B 


THEORY  DERIVATIONS 

1,  The  following  Is  the  derivation  for  equation  (1)  of  the  theory 
section  as  given  by  Stepanoff  /13/.  Newton's  second  law  is  as  follows: 

tr  —  d(AAV)  (B-  i 

dt: 

Where:  F  =  Force,  lbf. 

V  =  Velocity,  ft/sec. 
t^  =  Time,  sec. 

M  =  Mass,  slugs. 


* —  X-  <3  t  • 


The  velocity  through  the  pipe  (V  )  before  valve  closure  is  constant.  The 
force  (increase  in  pressure  at  the  valve)  is  equal  to  the  time  rate  change 
of  momentum  of  the  moving  mass.  The  pressure- force  relationship  is  as 
follows : 

F  =  A  v  H3  (b-2) 

At  time  t.  after  closure  the  portion  of  the  fluid  in  the  length  of  pipe,  X, 
will  have  its  velocity  reduced  to  zero,  then: 

F=A<TH,=  V.  <»-3> 

2 

Where:  A  =  Pipe  internal  cross  sectional  area,  ft 

Y  =  Specific  weight  of  fluid  flowing,  lbm/ft^t 
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H  =  Head  rise  above  initial  no  flow  static  head,  H  ,  ft  of 
y  o 

J  water 

g  =  Proportionality  constant  in  Newton's  second  law  of  motion, 
(lbm-ft)/ (lbf-sec2) 

X  =  Distance  to  any  position  on  the  pipe  measured  from  the 
valve,  ft. 

t^  =  Time  after  closure,  sec. 

Vq  =  Velocity  of  fluid  in  the  pipe  before  valve  closure,  ft/sec. 

The  term  A  —  in  equation  (B-3)  is  the  mass  flow  rate.  In  equation 
“o 

(B-3),  X/ 1 ^  is  the  velocity  of  the  pressure  pulse,  a.  Substituting  a  and 
cancelling  the  Y  A  term  on  both  sides  of  the  equation  (B-3)  gives: 


Li  ^  Vo 

=  a  -57  (1) 

2.  The  derivation  of  equation  (2)  in  the  theory  section  as  given  by 
Stepanoff  / 1 3/ ,  is  as  follows.  After  instantaneous  valve  closure  a  region 
of  compressed  water  is  created  adjacent  to  the  valve  and  extends  a  distance 
X  at  time  t^  after  closure.  The  pipe  in  this  region  expands  under  pres¬ 
sure  from  D  to  D- .  From  the  instant  of  closure  until  t.  a  volume  of  wa- 
o  1  i 

ter  equal  to  AV^t^  passes  through  the  original  cross  section  at  X.  There 
is  no  outlet  for  the  water;  therefore,  the  volume  must  be  accommodated  by 
expanding  the  pipe  radially  and  by  compressing  the  water. 


hocus 
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The  increase  in  volume  caused  by  the  compression  of  the  water  in  the  pipe 
is  obtained  from  the  definition  of  bulk  modulus  as  follows: 


dP 


Kw  cHG 


-  Gw 


w 


p_ 

AGW 


(B-4) 


Where:  K  =  Water  bulk  modulus,  lbf/ft  . 

w 

3 

=  Volume,  ft 

AGy  =  Incremental  change  in  volume,  ft' 


P  =  Incremental  change  in  pressure  caused  by  the  velocity  head 
changing  to  pressure  head,  lbf/ft  • 


Solving  equation  (B-4)  for  AGW  gives: 

=  ~G*  £- 
yv 


(B-5) 


but  the  original  volume  G„  is  equal  to  AX.  Substituting  AX  for  G„  gives: 

w  W 


AGw  =  -AXJL 

*  w 


(B-6) 


The  negative  sign  indicates  that  AG  decreases  with  increasing  pressure. 
The  decreased  volume  is  used  to  accommodate  the  additional  mass  that 
enters  the  section  of  Pipe  X  feet  long.  The  volume  made  available  by  the 
compression  of  the  original  water  is: 

 A  X  P 


AGW  - 


K 


(B-7) 


w 


The  increased  volume  caused  by  the  expanding  of  the  pipe  is: 


A  Gp  =  f(P.*-D;)X 


(B-8) 


By  equating  the  additional  volume  of  water  that  must  be  accommodated 
(AV^t^)  to  the  volume  that  the  pipe  expands  plus  the  volume  that  the 
water  compresses  the  following  is  obtained: 
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Av.ti  -  -45  +  ^(o;--o/jx 


For  the  pipe 


POe 


P  D0  =  2Se 


Where: 

e  =  Pipe  wall  thickness,  ft, 
Hooke's  law  for  the  pipe  expansion  is: 

D.-O, 


-E( 


Dc 


Where:  E  =  Modulus  of  elasticity,  lbf/ft  . 

Combining  equations  (B-10)  and  (B-ll)  gives: 


(B-9) 


(B-10) 


(B-ll) 


P  Do 

2  aE 


P,-Oo 

Do 


(  B  - 1 2  ) 


2  2 

In  equation  (B-9)  -D  may  be  written  (D^+Dq) (D^-Dq) . 

Since  Dj'D0  is  very  small  compared  to  or  Dq,  the  term  (D^+D^)  may  be 

written  2Dq.  Then  the  last  term  of  equation  (B-9)  becomes: 

_  ~x.fr  P XP.-0-)  (B.13> 

4  ~  z 

Multiplying  both  sides  of  equation  (B-12)  by  gives: 


d(-o.  =  JLSfL 

2  e  E 


(B  - 1 4 ) 


Substituting  equation  (B-14)  into  equation  (B-13)  and  substituting  A  for 

gives 

x  it  ( p Dc1)  _  2L  a  PqP 
4-  eE 


the  cross  sectional  area  term, 


(B- 15) 


Substituting  equation  (B-15)  for  the  last  term  in  equation  (B-9)  gives: 
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V„  A  t; 

V.t  i 


X  PA 

xp(± 


H- 


x  Oo  pa 

e.  E 

Po  \ 

eEJ 


(B-16) 


The  term  X/t^  is  the  wave  velocity,  a.  Substituting  and  solving  for  a 


gives : 


Q  = 


Vo 


(  -L  4-  JBi) 
V  Kw  eE/ 


(B- 17) 


If  equation  (1)  is  multiplied  by  the  specific  weight  of  water  (  "$"  ) ,  the 


following  is  obtained: 


1 1  v  _  o _  ^  v» 

2TW  —  r  —  cj0 


(B-18) 


Solving  (B-18)  for  Vq  gives: 


V.  = 


Pc 


Q  r 


w 


Substituting  (B- 19)  into  (B -17)  and  letting  D  equal  and  J ^  equal 
gives : 


(B-19) 

X, 


a  = 


+•  ek) 


(2) 


The  maximum  theoretical  pressure  rise,  (P^p)  *s  given  by  substituting  equa¬ 


tion  (2)  into  equation  (B-18). 

P.X  = 


V.  V, 


tTW) 


(A) 


The  experimental  results  were  expressed  as  P/P^, 


P  _  P  '/p-C^L  *  ~fe) 


Vo 


(B-20) 


The  term  P  is  the  experimental  pressure  rise  above  the  no  flow  static  pres¬ 
sure  * 
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3.  The  following  is  the  derivation  of  equations  (3),  (3a),  (5),  (5a), 

(6),  and  (6a)  of  the  theory  section.  These  derivations  are  based  upon 
work  done  by  Karplus  /4/.  His  equations  were  modified  to  make  them  apply 
to  the  water  hammer  model.  The  following  is  a  form  of  equations  (2)  for 
an  air-water  mixture. 


(B-21) 


Where:  Subscript  c  represents  the  air-water  mixture. 


c 


2 

mixture  compressibility,  ft  /lbf. 


The  density  of  the  mixture  ( 


c 


)  is  given  by: 


(B-22) 


3 

Where:  G  =  Volume^ ft  . 


M  =  Mass,  slugs. 


subscripts 


a  =  Air 


W  =  Water 


Because  mass  is  the  product  of  density  and  volume ,0  becomes : 


(B-23) 


The  volume  fraction  of  air  is  defined  as: 


(B-24) 


The  volume  fraction  water  is  defined  as: 


^  =  -§? 


(B-25) 


Substituting  (B-24)  and  (B-25)  into  (B-23)  gives; 


(B-26) 
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Compressibility  is  defined  as  follows: 


X 


c 


c 

—  (d  Ga  +•  olGw) 

GcdP 

Gt  dP 

olGw  (Gw \ 

**  / ^ 

G  w  d  P  VGc  / 

<*°l  P  VG«/ 

Substituting  equations  (B-24)  and  (B-25)  into  (B-27)  gives: 


X= 


(B-27) 


(B-28) 


Substituting  equations  (B-26)  and  (B-28)  into  equation  (B-21)  gives: 


Q- *=[(!-*)/„+  x/>„][xX,  +0-x)Xw  + 


(B-29) 


Rearranging  equation  (B-29)  gives: 


For  isothermal  conditions  and  assuming  the  ideal  gas  law  applies 


% 


IQ 


ASS 


Where:  ^ABS  =  Absolute  Pressure  of  system,  lbf/ft  . 

2 

^  =  Isothermal  compressibility  of  air,  ft  /lbf 

The  velocity  of  sound  in  air  is  given  as: 

Where  %a=  Adiabatlc  compressibility  of  aii^  ft  /lbf. 

For  adiabatic  conditions  and  assuming  the  idea'l  gas  law  applies: 

^  —  1 
A  Aft  A>  i 


A  R 


A8S 


(B-31) 


(B-32) 


(B-33) 
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Where : 


k  =  ratio  of  specific  heat  of  air  at  constant  pressure  to  that 


at  constant  volume,  1.403  cor  air 

By  combining  equations  (B-31)  and  (B-33),  the  following  is  obtained: 

<»-34> 

The  velocity  of  sound  in  water  is  given  by: 

=  /w  “K*, 

For  isothermal  conditions  substitute  (B-32),  (B-33),  (B-34),  and  (B-35)  into 
(B-30).  This  gives: 


Q 


-2 
X  ~ 


+■ 


(3) 


Using  the  adiabatic  conditions  for  equations  (B-30)  gives  the  following 
results : 


(3a) 


For  isothermal  conditions  for  an  air-water  mixture  equation  (B-18)  becomes: 


Vo  Qr 


(B-36) 


Where:  a^  is  given  by  equation  (3). 

Substituting  equation  (B-36)  into  (B-20)  gives: 


(6) 


Where :  a^.  is  given  by  equation  (3). 
For  adiabatic  conditions: 
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(6a) 


Where:  a  is  given  by  equation  (3a). 

A 
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APPENDIX  C 


EQUIPMENT 

1.  The  head  tank  consisted  of  a  tank  which  was  two  feet  in  diameter  at  the 
lower  end  and  one  foot  ten  and  one-half  inches  at  the  upper  end.  The  upper 
part  of  the  tank  was  light  sheet  metal  and  the  lower  part  of  the  tank  was 
plate  steel.  Two  inch  overflow  lines  were  provided  at  approximately  two 
foot  intervals  as  is  shown  in  Appendix  A.  These  overflow  lines  were  stan¬ 
dard  two  inch  galvanized  steel  pipe.  Each  overboard  line  had  a  two  inch 
gate  valve  in  the  line.  All  the  overboard  lines  ran  into  a  common  line  which 
discharged  into  a  collecting  tank  open  to  the  atmosphere.  This  collecting 
tank  had  a  drain  which  discharged  into  the  main  sewage  system,  A  Weston 
thermometer  located  near  the  bottom  of  the  tank,  was  used  to  measure  the 
inlet  water  temperature. 

The  main  run  of  pipe  consisted  of  two  inch  extra  strong  black  steel 

pipe.  The  out  side  diameter  was  2,375  inches,  the  inside  diameter  was  1.939 

inches.  The  wall  thickness  was  0.218  inches.  Three  sections  of  pipe  were 

welded  together.  Care  was  taken  to  insure  good  butting  of  the  individual 

♦  * 

sections.  At  the  transducer  locations,  transducer  fittings  one  inch  in 
diameter  and  one  inch  high  were  welded  into  the  pipe.  These  fittings  had 
standard  14  mm  by  1,25  automotive  spark  plug  threads. 

The  water  we*gh  system  consisted  of  a  Toledo  Scales  Company  platform 
scales  with  a  capacity  of  1600  pounds  and  a  steel  plate  rectangular  tank, 

20  inches  by  28  inches  by  20  inches.  A  3  inch  by  3  inch  section  was  cut 
out  of  the  top  to  permit  entrance  of  the  discharge  pipe.  This  tank  had 
two,  2  inch  gate  valves  for  discharge.  They  discharged  into  drain  lines  to 
the  sewage  system. 
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The  following  scale  information  is  provided. 


Manufacturer 


Toledo  Scales  Company 
Toledo,  Ohio 


Model 


31-1811FC 


Serial  Number 


2180-0-225 


Capacity 


0-1600  pounds 


An  electric  timer,  which  measured  to  the  nearest  one  tenth  cf  a 
second,  was  used  to  measure  the  time  to  fill  the  weigh  tank  with  water. 

2.  The  air  injection  system  is  shown  schematically  as  part  of  Fig.  1. 

A  photograph  of  the  air  injection  system  is  shown  in  Fig.  2.  The  source  of 
air  was  from  the  125  psi  air  system  permanently  installed  in  Building  500, 
at  the  United  States  Naval  Postgraduate  School.  This  system  was  reduced 
to  approximately  fifteen  psi  by  a  regulating  valve.  One-quarter  inch  vinyl 
tubing  was  used  from  the  regulating  valve  to  a  pipe  manifold.  A  copper- 
constantan  thermocouple  was  installed  in  the  pipe  manifold  to  measure  the 
temperature  of  the  air.  A  Rubicon  Company  catalog  number  2732  potentiometer 
was  used  to  measure  the  thermocouple  output.  This  is  item  1  in  Fig.  2. 

To  measure  the  pressure  of  the  air  a  50  inch  mercury  manometer  was  used. 

A  one-quarter  inch  vinyl  tubing  line  was  tapped  into  the  air  manifold.  This 
line  was  connected  to  the  mercury  manometer.  The  other  side  of  the  manometer 
was  open  to  the  atmosphere.  Item  4  of  Fig.  2  shows  this  manometer.  From 
the  manifold  a  one-quarter  inch  vinyl  line  led  to  the  rotometer.  The  follow¬ 
ing  data  is  provided. 

Manufacturer  Fisher-Porter  Company 

Flowrotor  tube  number  2-F-1/4-20-5/70 

Float  glass 


Range  (14.7  psi,  70°F) 


0-0.325  ft3/min. 
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The  manufacturer's  calibration  curve,  converted  from  cubic  centimeters  per 
minute  to  cubic  feet  per  minute,  was  used.  Fig.  14  is  a  copy  of  this 
curve.  Four  points  on  this  curve  were  checked  using  American  Meter  Com¬ 
pany  wet  gas  meter  Al-20-11472. 

The  flow  rate  was  corrected  to  14.7  psia  and  70°R  for  each  run  using 
the  following  relationship  suggested  by  the  manufacturer  / 2/ . 

Q  =  ft'Vmin  at  70°F,  14.7  psia 

3 

=  pressure  at  which  measurement  was  made,  psia 

=  temperature  at  which  measurement  was  made,  °R. 

3 

Q  =  flow  rate  in  ft  /min  taken  from  the  curve  of  Fig.  14. 

This  curve  was  entered  with  the  rotometer  reading. 

Item  number  2  in  Fig.  2  shows  the  rotometer.  A  one-quarter  inch  vinyl 
tubing  connected  the  rotometer  with  a  manual  control  valve.  The  control 
valve  was  used  to  control  the  air  flow  rate.  From  the  control  valve  a 
one-quarter  inch  vinyl  tube  led  to  ;  manifold  section  in  the  pipe.  A 
stop  valve  was  provided  for  each  of  these  separate  lines.  Item  5  of  Fig. 

j 

2  shows  the  equipment  used  for  the  injection  of  the  air  into  the  pipe. 

Three  separate  rings  of  holes  were  used  to  inject  the  air.  One  ring  con¬ 
sisted  of  two  holes,  each  with  a  diameter  of  0.04  inches.  The  second  ring 
consisted  of  four  holes,  each  with  a  diameter  of  0.07  inches.  The  third 
ring  consisted  of  eight  holes  with  diameters  of  0.07  inches.  These  holes 
were  equally  spaced  around  the  circumference  of  the  pipe.  A  plexiglass 
ring,  four  inches  in  diameter,  with  a  circumferential  groove  for  air  injec¬ 
tion  was  provided  for  each  ring  of  holes.  These  rings  were  sealed  with  "0" 
rings.  Fig.  5  shows  a  cross  sectional  view  of  this  system.  The  system  was 


50 


designed  so  that  any  one  of  the  rings  could  be  used.  However,  it  was 
not  necessary  to  do  this.  All  experimental  work  was  done  with  air  enter¬ 
ing  all  three  rings  of  holes. 

3.  The  valve  and  closure  system  is  shown  in  Fig.  3.  The  quick  closing 
valve  was  a  modified  standard  Crane  two  inch  quick  closing  gate  valve. 

The  valve  was  modified  by  screwing  on  a  steel  disc  two  inches  in  diameter 
and  one  inch  high  on  the  top  of  the  stem.  This  disc  was  the  point  of  im¬ 
pact  for  the  falling  weight  which  closed  the  valve.  A  mechanism  was  de¬ 
signed  to  hold  the  valve  closed  when  it  was  in  the  shut  position.  A  stand 
that  was  completely  isolated  from  the  piping  system  was  used  to  stop  the 
valve.  By  using  this  method,  the  full  impact  of  the  stopping  valve  was 
not  transmitted  to  the  piping  system.  A  seven  pound  cylindrical  steel 
bar  one  inch  in  diameter  by  24  inchs  long  was  used  as  the  weight  for  clos¬ 
ing  the  valve.  This  weight  was  dropped  through  an  aluminum  pipe,  one 

and  one-quarter  inches  in  diameter.  This  aluminum  pipe  contained  numer¬ 
ous  air  holes  to  allow  the  weight  to  fall  freely.  Appendix  G  contains  pre¬ 

liminary  calculations  used  to  estimate  the  closure  time.  The  weight  was 
dropped  from  nine  feet  for  all  runs.  The  closure  time  of  the  valve  was 
measured  with  an  oscilloscope.  The  time  sweep  of  the  oscilloscope  was 
triggered  by  a  photoelectric  cell.  A  contact  completed  a  45  volt  direct 
current  circuit  when  the  valve  was  closed.  This  voltage  was  put  across 
the  vertical  plates  of  the  oscilloscope.  A  Polaroid  time  exposure  was 
made  of  the  oscilloscope,  the  closure  time  was  taken  from  the  photograph. 
Polaroid  3000  speed  film  was  used. 

The  quick  release  mechanism  consisted  of  a  steel  pin  which  held  the 
weight  nine  feet  above  the  valve.  A  line  was  attached  to  this  pin.  When 
it  was  desired  to  drop  the  weight,  the  pin  was  pulled  and  the  weight  fell. 
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A  nylon  string  was  attached  to  the  weight.  This  was  used  with  a 


pulley  to  pull  the  weight  to  the  desired  height  above  the  valve. 

4.  Selected  characteristics  of  the  pressure  sensing  equipment  are  given 
below.  For  complete  details  refer  to  manufacturer’s  instruction  book  / 5/ . 

Kistler  401  Pressure  Transducer 


pressure  range 

0.0-3000  psi 

overload  capacity 

507. 

linearity 

17. 

natural  frequency 

48 s 000  cps 

maximum  continuous 
temperature 

500°F 

maximum  intermittent 
temperature 

3000 °F 

length 

2.04  in. 

diameter 

0.625  in. 

weight 

1.75  oz. 

serial  number 

#1044 

Kistler  Model  651  Piezo  Calibrator 


input  voltage ,  maximum 

0.6  volts 

output  voltage,  maximum 

0. 3  volts 

voltage  gain 

0.5 

input  impedance 

,„14 

10  ohms 

output  impedance 

lO^ohms 

frequency  range 

0-10  kc 

linearity 

17. 

dimensions 

7  x  10  x  6  in. 

weight 

4.7  lbm 
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Hewlett  Packard  Model  130  Oscilloscope 


internal  sweep  accuracy  5% 

Vertical  and  Horizontal  Amplifiers 
frequency  response  0-300  kc 

sensitivity  lmv/cm-20  volts/cm 

accuracy  570 

input  impedance  1  megohm 

Cathode  Ray  Tube 

power  requirements  115/230  volts,  AC 

50/400  cps,  approximately 
175  watts 

dimensions  9  3/4  in.  by  15  in.  by  21-1/4  in. 

weight  39  lbm 

Kistler  transducer  model  401  serial  number  1044  was  calibrated  using 
an  Ashcroft  dead-weight  tester*  The  system  used  for  calibration  was  the 
same  as  shown  in  Fig.  4  except  the  transducer  was  inserted  into  the  14  mm 
dead-weight  adapter  fitting  instead  of  the  pipe.  Sweeney  / 1 4/  describes 
the  procedure  for  using  the  dead-weight  tester.  The  important  thing  in 
calibrating  was  to  rotate  the  weights-  when  loading  the  transducer.  The 
transducer  was  rapidly  loaded  so  that  there  would  be  no  error  due  to 
charge  leakage  from  the  crystal.  The  transducer  was  calibrated  for  the 
calibrator  range  scales  of  10  and  lOf.The  10  range  covers  0-125  psig  and 
100  range  covers  0-1250  psig.  Fig.  15  and  16  show  these  calibration 
curves.  These  curves  show  the  readings  on  the  calibrator  for  a  given 
pressure.  It  was  necessary  to  determine  the  number  of  calitration  divisions 
that  correspond  to  a  centimeter  deflection  on  the  oscilloscope.  Knowing 
this  and  using  the  calibration  curves*  the  psig/cent imeter  on  the  oscillo¬ 
scope  was  calculated.  For  the  oscilloscope  used  for  this  experimental 
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PIG.  15*  Transducer  Calibration  Curve,  0  -  125  ps’ig 
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PIG#  16#  Transducer  Calibration  Curve,  Q-1250  psig 


work,  the  following  values  were  obtained  and  used.  These  values  are  for 


Kistler  model  401  serial  number  1044  pressure  transducer,, 


Calibrated 

range 

Osci llcseopa 
scale 

Calibration 

divisions/ 

10  cm  oscillo¬ 
scope  deflection 

psig/ cm 

10 

5  mv/ cm 

19.7 

3.0 

10 

10  mv/cm 

39.1 

6.0 

10 

20  mv/cm 

77.2 

11.8 

100 

5  mv/cm 

17.2 

26.6 

100 

10  mv/cm 

34.4 

54.4 

100 

20  mv/cm 

67.8 

106.9 
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APPENDIX  D 


DETAILED  PROCEDURE  FOR  A  RUN 

The  detailed  step  by  step  procedure  for  each  run  was  as  follows: 

1.  Record  atmospheric  temperature  and  pressure. 

2.  Fill  the  tank  to  the  level  desired. 

3.  Open  quick  closing  stop  valve. 

4.  Adjust  water  inlet  to  the  tank  until  the  overflow  at  the 
desired  level  is  just  a  trickle. 

5.  Adjust  the  air  rotometer  to  desired  flow  rate  with  the 
^control  valve  immediately  following  the  rotometer  in  the  air 

inlet  system. 

6.  Insure  weigh  tank  discharge  valves  are  completely  open  and 
that  the  tank  is  as  low  as  possible. 

7.  Set  electric  timer  to  zero. 

8.  Close  weigh  tank  discharge  valves. 

9.  Start  timer  and  record  initial  weight  on  scales  when  the 
timer  is  started. 

10.  Fill  weigh  tank  until  it  is  nearly  full. 

11.  Record  weight  and  time  at  nearly  full  position. 

12.  Open  weigh  tank  discharge  valves. 

13.  Lift  weight  to  nine  foot  position  and  insert  quick  release 
pin. 

14.  Record  inlet  and  outlet  water  temperatures, 

15.  Record  rotometer  setting  (read  ball  top). 

16.  Record  air  mercury  manometer  readings. 

17.  Record  air  temperature. 
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18.  Set  calibrator  range  to  desired  position.  Insure 

the  meter  is  adjusted  to  1,0.  Allow  a  ten  minute  warm  up 
period  for  the  calibrator  and  the  oscilloscope, 

19.  Set  oscilloscope  on  desired  vertical  scale. 

20.  Adjust  oscilloscope  sweep  rate  to  the  desired  time  scale. 

21.  Insure  the  SYNC  switch  is  set  to  external  trigger  (EXT)  on 
the  oscilloscope. 

22.  Insure  that  battery  is  properly  connected  to  trigger 
circuit  of  the  oscilloscope. 

23.  Set  the  horizontal  sensitivity  of  the  oscilloscope  to 
internal  sweep. 

24.  Check  system  by  checking  the  calibration  reading  for  ten 
centimeters  deflection  on  the  vertical  scale  of  the  oscillo¬ 
scope. 

25.  Adjust  vertical  and  horizontal  position  until  the  oscillo¬ 
scope  trace  is  in  desired  position. 

26.  Set  camera  on  time  exposure. 

t  ¥ 

27.  Open  camera  shutter. 

28.  Close  valve  by  dropping  the  weight  with  the  quick  release 
mechanism. 

29.  Close  camera  shutter. 

30.  Open  quick  closing  valve. 

31.  Develop  Polaroid  photograph. 
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HORIZ  NTAl.  EXPERIMENTAL  OATA  SHEET 


APPENDIX  F 


SAMPLE  CALCULATIONS 

Sample  calculations  have  been  included  for  run  number  3-9-10.  The 
data  for  this  run  has  been  included  on  the  sample  data  sheet  in  Appendix 
E. 

PATM  =  30 • 264  in*  HS 

t«™  =  76. 0°F 
ATM 

c  =  -0.129  Barometer  temperature  correction 
Calculate  atmospheric  pressure  in  psia 

(pATM(in  «8)  +  c) (0.4913)  =  pATM(Psia) 

(30. 264-0. 129)(0. 4913)  =  14.81  psia 

where  0.4913  is  the  conversion  factor  for  converting  pressure  in  inches 
of  Hg  to  psia. 

Calculate  absolute  pressure  of  air. 


(Pl3  +  ?2a)  (0'4913)  +  PATM  =  pa(Psia> 

P^a  and  are  the  inches  of  mercury  measured  on  each  leg  of  the  air 

manometer. 

Pla  =  12’62  in*  Hg 
P2a  =  12.55  in.  Hg 

PATM  *  U'81  pSi“ 

=  pressure  at  which  air  was  measured,  psia 
(12.62  +  12.55)(0.4913)  4-  14.81  =  27.18  psia 
Calculate  air  flow  rate  at  14.7  psia  and  70°F. 
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where: 


Q  =  air  flow  rate  at  14.7  psia  and  70°F,  ft  /min. 

d 

Q  =  uncorrected  air  flow  rate  taken  from  Fig,  14, 
u 

Enter  Fig,  14  with  rotometer  setting. 

=  Absolute  temperature  at  which  air  flow  rate  was  measuredj 


°R  / - 

/  Sio  \ 

Qa  --  °-325  J  LtjzJ 


=  0.442  ft^/min 


Calculate  mass  flow  rate  of  air 


™a  =  Qa  X* 


where : 

=  Mass  flow  rate  of  air,  lbm/min 

y  =  Specific  gravity  of  air  at  14.7  psia  and  70°F. 
a 

ma  =  (0.442)  (0.0752)  =  0.03403  lbm/min 
Calculate  volume  flow  rate  at  inlet  water  conditions 


=  Qa(Jr£)  (yk) 


where 


Q  =  Air  volume  flow  rate  at  conditions  of  water  inlet 
aw 

=  Absolute  pressure  of  water  at  air  inlet,  psia 
=  Absolute  temperature  of  air  at  water  inlet,  °R 


Q  =  0.442 
aw 


(J±l  )(. 

V  18.245/ L 


£i±  ) 

S3oJ 


=  0.3463  ft  /min 


Calculate  the  water  mass  flow  rate 

Yrt  .=  foo) 

v-J  4.  • 

• 

Yfl =  Water  mass  flow  rate,  lbm/min 

•  '  W 

VY1W  =  Weight  of  water  weighed  in  time  t^,  lbm 
t^  =  time  to  measure  pounds  of  water,  sec 
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mtT  =  UlPHpp)  =  534  lbm/min 

w  21.4 

Calculate  water  volume  flow  rate 


where : 


*W 


ra  w 


3 

Q._  =  Water  volume  flow  rate,  ft  /min 
W 

Yw=  Specific  weight  of  water,  lbm/ft 

QW  =  ^624“  =  8-56  ft3/min 

Calculate  the  volume  percent  air  in  mixture 


Q 

Q 


aH<100> 

aU+% 


where : 


=  Volume  percent  air  at  water  inlet  conditions  % 

"  ’  (0.3463  +  5.56)  ~ 

Calculate  the  weight  percent  air  in  the  mixture 

vm) 

Wl,^  +■ 


where : 


n 


WT 


=  Weight  percent  air  in  air  water  mixture,  % 


tf  (0.03403) (100) 

UW7  ~  (0703403  +  534)  =  0*006377- 


Calculate  the  maximum  pressure 
P  =  d#C-J 

P  =  Pressure  above  the  static  no  flow  pressure,  psi 
C  =  Psi/cm  of  oscilloscope  deflection,  psi/cm 
d  =  Deflection  of  oscilloscope,  cm 

J  =  Correction  for  velocity  head  at  the  location  when  the 
oscilloscope  scale  was  zeroed,  psi. 
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To  determine  the  J  correction,  the  difference  in  static  pressure  with  flow 
and  with  no  flow  was  experimentally  measured  for  each  head  at  each  loca¬ 
tion.  It  was  necessary  to  make  this  correction  because  the  oscilloscope 
was  zeroed  for  each  run  when  there  was  velocity  head  in  the  pipe. 

P  =  (3. 1) (6.0)-2, 1  =  16.5  psi 
Calculate  pressure  ratio  P/P^ 


P/P 


NT 


P/P 


/x  -k,\i  i  ) 

.,/£7I7K)  ’• 


NT  w  v 


where : 


Vo*w 


Vq  =  Velocity  of  the  water  through  the  pipe  with  no  air, 

determined  experimentally,  ft/sec 

2 

=  Density  of  the  water,  slugs/ft 

2 

Bulk  modulus  of  water,  lb/ft 
D  =  Inside  diameter  of  pipe,  in. 
e  =  Pipe  wall  thickness,  in. 

2 

E  =  Modulus  of  elasticity  of  pipe  lbf/ft 

=  Theoretical  maximum  head  rise  above  the  initial  static 
head,  equation  l,ft  of  water 


) 


1.51$ 


p/p  =  UdhL.  (3l0'°oOXl44,)  10^44) 


P/P„^  =  0.0445 
NT 
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APPENDIX  G 


VALVE  CLOSURE 

The  water  hammer  model  of  Fig.  1  was  used  to  generate  the  pressure 
pulses.  For  the  maximum  pressure  to  be  theoretically  possible,  the  model 
required  a  valve  which  would  close  in  less  than  a  LkzA 

seconds. 

9 

where:  L  =  The  length  of  the  pipe  from  the  tank  to  the  valve,  ft. 

X  =  The  distance  from  a  given  location  to  the  valve,  ft. 

a  =  The  velocity  of  the  pressure  pulse  in  the  pipe,  ft/sec. 

The  speed  of  the  pressure  pulse  in  pure  water  in  a  two  inch  extra  strong 

steel  pipe  is  4580  ft/sec.  For  the  maximum  pressure  as  given  by  equation 
(1)  to  be  possible  at  transducer  location  two,  the  valve  had  to  close  to 
less  than  0.0128  seconds.  If  there  was  0.05  percent  by  volume  air  at  the 
water  inlet  conditions  tie  closure  time  required  at  transducer  location 
two  for  the  maximum  pressure  to  be  possible  was  0.0351  seconds.  With  ten 
percent  by  volume  air,  the  closure  time  required  at  transducer  location 
two  was  0.438  seconds  or  less.  Since  the  pressure  wave  travels  fastest 
in  water  with  no  air,  this  was  the  limiting  condition.  Therefore,  to  use 
transducer  location  two,  a  closure  time  of  less  than  0,0128  seconds  was 
required. 

A  free  falling  seven  pound  weight  starting  with  zero  initial  velocity 
and  falling  nine  feet  would  have  a  velocity  calculated  as  follows: 

V  =  V  2gs 

V  =  v/2 (32 , 2) (2)  =  24.0  ft/sec 

2 

where  g  =  acceleration  of  gravity,  ft/sec 
s  =  distance  the  weight  falls,  ft 
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V  =  velocity  of  the  falling  weight,  ft/sec 
With  a  24.0  ft/sec.  velocity,  the  time  of  closure  would  be: 


where  t^  =  time  for  valve  gate  to  travel  Sv  distance,  sec 
Sv  =  distance  valve  gate  travels  for  closure,  ft 
V  =  velocity  of  valve  gate 
tt  =  (1/ 6) (1/24)  =  0.0069  sec. 

Assuming  an  plastic  collision,  the  following  calculations  are  made. 

M1V1  =  (M2  +  M1)V2 
M1 

V  =  — - -  V 

2  M1+M2  1 

Assuming  M2  =  1.75  lbm  =  1.75/32.2  slugs 

Mj  =  7  lbm  =  7/32.2  Slugs 

/_7 _ 1  (24)  =  19.2  ft/sec. 

V2  =\(7+l. 75)^ 

The  closure  time  for  the  valve  would  be: 

fci  =  Si/V2  =  (6) (19.2)  =  °*0087  sec* 


The  packing  on  the  valve  was  loose  enough  that  the  valve  closed  freely 
with  very  little  resistance.  The  weight  of  the  valve  stem  and  gate  was 
1.75  lbm. 

Experimentally  the  closure  time  was  measured  as  0.007  seconds  with  an 

uncertainty  of  0.003  seconds.  The  method  of  measuring  the  closure  time 

♦ 

lias  been  discussed  in  Appendix  C. 
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APPENDIX  H 


UNCERTAINTY  ANALYSIS 

The  uncertainty  for  runs  3-9-2  and  3-9-10  are  included  here  as  re¬ 
presentative  examples  of  the  uncertainty  of  the  experimental  work.  Since 
the  magnitude  of  the  variables  differs  between  runs,  these  values  are  to 
be  considered  only  as  representative.  The  method  of  analysis  is  that  pro¬ 
posed  by  Kline  and  McClintock  / 6/  for  single  sample  data. 


Run  Number 

3-9-2 

% 

3-9-10 

7« 

Item,  Units 

Value  +  Uncertainty 

Value  +  Uncertainty 

m  ,  lb/min 

V 

563.6  +  13.6 

2.4 

534  +  12 

2.2 

3 

Q  ,  ft  /min 
w 

9.03  +  0.22 

2.4 

8.56  +  0.19 

2.2 

m  ,  lb 
w* 

170+  4 

2.4 

190  +  4 

2.1  •  . 

water 
t.,  sec 

18.1  +  0.1 

0.6 

21.4  +  0.1 

0.5 

T  ,  °R 
w 

516  +  1 

0.2 

518  +  1 

0.2 

P  ,  psia 
w 

18.26  +  0.07 

0.4 

18.26  +  0.07 

0.4 

3 

Q  ,  ft  /min 
u 

0.040  +  0.004 

1.0 

0.325  +  0.004 

0.1 

3 

Q  ,  ft  /min 

cl 

0.0562  +  0.006 

10.0 

0.442  +  0.0057 

1.3 

3 

Q  ,  ft  /min 
aw 

0.04403  +  0.00439 

9.9 

0.3463  +  0.0047 

1.4 

P,  psia 

51.9  ±  6 

10.6 

17.1  +  3 

17.5 

P/PNT' 

0.1121  +  0.0129 

10.6 

0.0368  +  0.0064 

17.5 

Patm>  "S 

30.264  +  0.002 

0.007 

30.264  +  0.002 

0.007 

Tatm’  °F 

76.0  +  0.5 

0.65 

76.0  +  0.5 

0.65 

ft,  * 

0.485  +  0.050 

10.2 

3.89  +  0.77 

0.8 
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